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Abstract 

We consider the problem of gravitational forces between point par- 
ticles on the branes in a Randall-Sundrum (R-S) two brane model 
with IZ2 symmetry. Matter is assumed to produce a perturbation 
to the R-S vacuum metric and all the 5D Einstein equations are solved 
to linearized order (for arbitrary matter on both branes). We show 
that while the gauge condition /ijs = 0,i = 0,1,2,3 can always be 
achieved without brane bending, the condition h^<^ = leads to large 
brane bending. The static potential arising from the zero modes and 
the corrections due to the Kaluza-Klein (KK) modes are calculated. 
Gravitational forces on the Planck (yi = 0) brane recover Newtonian 
physics with small KK corrections (in accord with other work). How- 
ever, forces on the TeV (7/2) brane due to particles on that brane are 
strongly distorted by large R-S exponentials. 



1 Introduction 



Over the past two decades, the possible existence of extra dimensions has 
played an important role in much theoretical analysis. Thus M-theory ex- 
ists in 10 or 11 dimensions where six are compactified, characteristically on 
a Calabi-Yau (C-Y) manifold. Examples of this are 10 dimensional (lOD) 
intersecting D-brane models (see e.g.^ and references therein), and 11 di- 
mensional Horava- Witten (H- W) M-theory [21 El IH E] ■ In the latter case the 
eleventh dimension is compactified with /Z2 symmetry so that one can 
think of space as bounded by the two lOD orbifold planes at ?/i = and 
2/2 = Trp with appropriate boundary conditions at yi and 2/2 to enforce the 
/Z2 symmetry. In H-W theory, only the gravity multiplet propagates in 
the bulk, while normal matter fields exist on the orbifold planes, the physi- 
cal plane is conventionally chosen to be at ?/i = 0, while supersymmetry is 
broken on the distant ?/2 plane Since vrp is (9(10) times larger than the 
Calabi-Yau radius (which is ^ l/Mg, where Mq ^ 3xlO^^GeV is the GUT 
mass) one can consider a 5D reduction of the theory with the effects of the 
C-Y space taken into account by moduli acting in the bulk, e.g. the C-Y 
volume moduli V = e'^ being the simplest example. 

The phenomenological 5D Randall- Sundrum (R-S) model [HI IZ| built on 
the same geometrical structure has received considerable analysis to see if 
such 5D models can reproduce standard 4D results. The initial question 
examined was whether such a model could reproduce the 4D Friedmann- 
Robertson- Walker (FRW) cosmologyjH iniliniiniliaiI31IIllliaiinilIIllIH|- 
In order to do this, an additional scalar field in the bulk was added (analogous 
to the field in H-W) to stabilize the vacuum metric with phenomeno- 
logically chosen potentials in the bulk and on the branes to fine tune the net 
cosmological constant on the brane to zero. Matter on the branes is then 
treated as a perturbation (as appropriate for the Hubble expansion era) and 
it was seen that indeed the 4D FRW cosmology could be recovered. A some- 
what surprising result was that the same result could be obtained for H-W 
theory only in the radiation dominated era |TH|(at least for the tree level the- 
ory). While here the potentials of H-W theory automatically guarantee the 
vanishing of the net cosmological constant without fine tuning, they are not 
free as in R-S models but are determined by the theory. One finds that the 
phenomenologically chosen R-S potentials allow the brane separation to be 
determined by the non-relativistic matter density pnr-, something that does 
not occur in H-W theoryjlSj. 
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A second question, which is the subject of this paper, is whether the 
5D models give rise to the correct 4D Newtonian gravitational potential 
(in leading order). To examine this, one considers point particles on the 
branes and calculates the gravitational forces between them. There is a large 
literature on this subject as well and on calculations of corrections arising 
from the presence of an additional dimension piTl 1 ^ 1 ^ 1 ^ 1 ^ 1 ^ 1 ^ 177] . 
This question brings up an additional feature of the R-S model. The vacuum 
metric of the R-S model has the form 

ds^ = e-^^^y\ijdx'dx^ + dy^ (1) 

where A(y) is an increasing function of y. (We use the notation i,j = 0,1,2,3; 
/i, z/ = 0,1,2,3,5 and = y.) Thus to account for the gauge hierarchy, the 
physical brane must be at ?/2 (where e"'^*-^^-' ^ 10^^^) while in the H-W theory 
one can assume that the physical brane is at ?/i = 0. Thus for R-S one needs 
to calculate the gravitational forces between the two particles at y2 (though it 
is also interesting to see what forces the theory predicts between one particle 
at i/i = and one at j/2 = vrp). Previous analyses have only examined forces 
between particles on the = brane^ and also neglect the effects of the 
Goldberger-Wise scalar stabilization field (j). For this case the function A(y) 
reduces to 

A{y) = P\y\ ; 2/1 - e < y < 2/2 - e ; e > (2) 

We examine here within this framework the general case of gravitational 
forces between particles on both branes as well as the size of the leading 
corrections due to the extra dimensions. We find that the force on a particle 
on ?/i = due to other particles on yi and ?/2 has a leading Newtonian form 
(in accord with previous work) though the Newton constant is different 
for the two cases. However, the Newtonian force on a particle at 2/2 due to 
another particle at ?/2 contains terms that grow exponentially with y2 which 
leads to an unsatisfactory theory. 

In carrying out these calculations, it is important to take careful account 
of "brane bending" effects. Thus we assume that matter is added on the 
branes ClS db perturbation to the vacuum metric 

ds"^ = e~'^'^y{r]ij + hij)dx'dx^ + hi^dydx' + (1 + h5)dy'^ (3) 

^R,ef.)22j considers a single brane model with the brane displaced from the origin. 
However this is different in that it does not have /Z2 boundary conditions imposed at 
y = yi and y = ya- 
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and then solve the Einstein equations to hnear order in /i^j,. The diffeo- 
morphisms of a 5D theory with / Z2 symmetry are those of R^xS^ which 
commute with Z2. This means that for the transformation 

x" x'" + e = x'' (4) 

one has that vanishes at the orbifold points, yi and 7/2: 

e(x\y,)^0^e(x\y2) (5) 

If one were to make a coordinate transformation with a non- vanishing 
then the branes become bent and this would create a complication when one 
imposes the Z2 boundary constraint on the branes, leading to the so-called 
brane bending effects. In previous analyses, the 5D Einstein equations were 
solved in Gaussian coordinates described by 

/i5^ = ; d^Kj = = ri'^Kj (6) 

In general, these cannot be achieved without brane bending occurring. We 
give here an alternate analysis which avoids these complications by making 
only coordinate transformations that satisfy Eq.(5). 

In Sec. 2 we give the metric decomposition along with our gauge choices 
which will allow us to solve the Einstein equations in the presence of matter 
on the branes without introducing brane bending effects. In Sees. 3 and 4 we 
explicitly solve the bulk Einstein equations to first order in the metric pertur- 
bation in the static limit, and subject these solutions to the brane boundary 
conditions. In Sec.4 we also find the poles of the transverse traceless piece of 
hij and show how the Kaluza-Klein modes produce corrections to the lead- 
ing static potential. In Sec. 5 we give the form of the Newtonian potential 
and show that Newton's constant differs depending on whether the gravita- 
tional force is due to particles on coincident or separate branes. Conclusions 
are given in Sec. 6 and the leading corrections to the Newtonian terms are 
discussed more fully in the Appendix. 
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2 Coordinate Conditions 



Following 4D analyses, it is convenient to decompose the metric into its 
transverse and longitudinal parts according to the ADM prescription^: 

hij = hjj^ + hjj + hij + hj^i (7) 

where hj^ is transverse and traceless^ {d'^hj^ = = h]) and hjj is transverse 
but (in general) possesses a trace: 

d% ^ ; {1^)1 ^f^O (8) 

We also decompose hi into transverse and longitudinal parts 

h^ = hj + ; d'hj ^ (9) 



and can write 



where 



hjj = l^ijf ; = Vij - Oij (10) 



o, - '4 (11) 

One can express each of the subparts in Eq.(7) in terms of hij. Thus 
taking the divergence of Eq.(7) gives 

d'hij = dp^h^ + D^/iJ (12) 

and 

d'd^hij = {aYh^ (13) 

Thus 

h% = O^fiuih'^' (14) 

and 

/if, = 0]Kk - OiPuh''' (15) 



^This decomposition was first introduced in |2H|- (Ref.|2ni is a more accessible recent 
reprint summarizing the ADM formalism.) The generalization to 4-space is trivial except 
for the ambiguity in defining l/D^ in Minkowski space. However, we will always be 
considering the static Newtonian limit here where V^, though the size (and correct 

definition of) the higher order dynamical effect are also of interest. 

■^Four dimensional indices are raised and lowered with the Lorentz metric ry^ . We use 
the notation hi j = djhi 
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Taking the trace of Eq.(7) and using Eq.(14) determines to be 

f = 7i''h,j (16) 

and since 

= hij - hjj - hij - hj,i (17) 
one has ^ 

hjj^ = ■Kik'Kjih'''' - -TTijTTklh''^ (18) 

Our general metric with matter on the branes has the form 

ds'^ = e-'^^^y^ {riij + hij)dx'dx^ + h^idydx' + (1 + h55)dy^ (19) 

In order to discuss clearly the issues of brane bending, we assume here that 
there exists a frame with no brane bending (as e.g. is required in H-W the- 
ory) i.e. in this frame the 4D branes are orthogonal to the fifth dimension. 
The vacuum metric of Eq.(l) is indeed a solution of the 5D Einstein equa- 
tions obeying the /Z2 boundary conditions in this frame (as can be seen 
below Eqs.(37, 38)). Since the perturbation due to matter, /i^j,, in Eq.(19) 
is general, we can assume there is a choice of h^^ that holds in a frame with 
no bending. We now ask what coordinate conditions can be imposed to sim- 
plify the metric but remain in a frame with unbent branes. We chose the 
coordinate condition 

/i5i = ; i = 0,l,2,3 (20) 

To verify that this can be achieved without any brane bending, we consider 
the infinitesimal gauge transformation 

h',,{x) = h,,{x) + ^5,^ + 6,5 - 2V%i^ (21) 

where F^^ are the Christoffel symbols for the vacuum metric {h^j,i, = 0). Hence 

h',^{x) = h^i{x) + es.i + e-^^e'^),, (22) 

Thus if in the initial gauge h^i were not zero, one can always set h'^^ — 
without brane bending e.g. by choosing ^5 = and solving for 6. Eqs.(20) 
are not a complete set of coordinate conditions (being only four in number), 
and one may ask what is the remaining gauge freedom that maintains 5h^i 
— 0. Prom Eq.(22), this implies 

= 6. + 6-^^(6^0,5 (23) 
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Writing = + where d^^f = 0, one finds 

= F^ix') (24) 
^5 = -('^''),5 (25) 

where is independent of y, and we have introduced the notation 

u,^e'^^y^U^\y) (26) 



It is understood in Eq.(25) that from Eq.(5) uj^{x^,yi) = = a;5(x*,y2) to 
avoid brane bending, which consequently also constrains {00^)^5. 

The remaining gauge freedom, Eqs.(24, 25), allows a residual gauge free- 
dom in and /irr: 



5hij = ujij + Uj^i - 2A'e ^^cus^- (27) 
(5/i55 = 2(e-2^a;5),5 (28) 

where A'{y) = dA/dy. Decomposing /i^ into its transverse and longitudinal 
parts one finds using Eqs.(24, 25) that 

5hl^ = ; 5hJ = oofix') (29) 
6f = 6A'e-2^(o;^),5 = -6A'e-^^u, (30) 
^(□^/i^) = 2n2cu^ + 2A'e-2^(o;^),5 (31) 

Since a;5(x\yQ,) = we see that is gauge invariant on the branes 

Sf{x\y^)=0 ■ a = 1,2 (32) 

Further, since u!5{x\y) is otherwise arbitrary, one may expand it around y = 

Ua ^ 

l^5{x\ y) = {y- VaW^ix', ya) + 2^y- yaf^l{x\ ya) + - (33) 

SO that 

mf{x\ y^)) = -QA'e-^^uj',{x\ y^) + ... (34) 
Hence one may choose (x*, ya) to set 

d,f{x\y^)=Q (35) 

which we will see below is a convenient further gauge choice. 
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3 Einstein Equations 



The action for our system is 

j d'xy/^gi—MlR + mlP^) (36) 
+ J2 [ <i''^^9{l^m^ + (-l)"6M|/3)5(y - y«) 

a=l,2 

where M5 is the 5D Planck mass, are the Lagrangians for point particles 
on the branes yi = and = Trp, and we have fine tuned the bulk and 
brane cosmological constants so that the net cosmological constant is zero. 
The vacuum equations of motion for the metric of Eq.(l) read 

^A" = {-ir^'f35{y-y^) (37) 
= (38) 

so that A = f3y for yi < y < y2 with / Z2 boundary conditions at the 
orbifold fixed points. The linearized first order equations read 

4'^^ = -Ji^ - Ig^uTiyaMy - ya) (39) 

5 a 

where T = g'^^T^^, and Tni,{ya) are the 4D stress tensors for Hence T^i 
= = T55. We will consider here only the static gravitational forces and so 
only Too need to be taken non-zero^. 
The 15 Einstein equations are then 

4P = ; i = 0,l,2,3 (40) 



^5 = -^E^'^^oo(y.)<^(y-y.) (41) 



and 



^ = -^1 E(^^^(^«) - hje-'^my - y^) (42) 



^In obtaining Eq.(39), is the first order part of R^^ omitting terms proportional 
to the eosmological constant, since these terms are precisely canceled by the cosmological 
constant sources on the right hand side as a consequence of the zero'th order equations 
Eqs.(37, 38). 
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The 10 equations Eq.(42) can be decomposed as 

v''R^^ = -^,Y.^oo{ya)S{y-ya) (43) 



Ml 



and 



Ml 



Eqs.(43) and (44) together pick out the "T" components and "L" component 
of Rij while Eq.(45) picks out the TT components. In this section we solve 
Eqs.(40, 41, 43, 44). Eq.(45) is discussed in Sec.4 below. 
Eq.(40) reads 

Rf^ = ^v^'d^idihki - diKk) + \A'd,h^, = (46) 

where A! = dyA. In terms of the decomposition of Eqs.(7) and (9), Eq.(46) 
reduces to 

- - 3A'dih,5 + a'^d^hj = (47) 

where = — 9q. In the static approximation, — > V^, the T part of 
Eq.(47) reads 

d^hj = (48) 

which says that hj is independent of y. We can thus use the gauge freedom 
of Eq.(29) to set hJ to zero 

hJ = (49) 
The remaining longitudinal part of Eq. (47) yields 

h5 = -^d,f (50) 

Note from Eqs.(28) and (30), the combination SA'h^^ + d^f^ is gauge in- 
variant. However in the gauge of Eq.(35), one has that h^^ vanishes on the 
branes, i.e. 

h55{x\ya)=0 (51) 
though in general it is non-zero off the branes. 
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Eqs.(41) and (43) allow us to determine and h^. One has for i?^5^ 

4? = {\dl - A'd,)ri'^hij + ^e'^D'h, + 2A'd,h, (52) 
and so Eq.(41) becomes, using Eq.(50) 

{ldl-A%){n^h^ - If) - "^o'd,f = --I3 J]e^^Too5(y -y.) (53) 

5 a 

In arriving at Eq.(53) we have made use of the fact that 
A"d,f ~ S{y - = in the gauge of Eq.(35). 

The full R^]'^ is 

e'^R^^ - {^dl-2A'd,)h,j-^A%jd,{7^^'h,r^k) (54) 
+ Tjijh^iA" - AA'^) + ^A'r],,d5h55 



Hence Eq.(43) becomes 



L_Y.e'^T,,{y^)5{y-y^) (55) 



3M| 



Subtracting Eq.(55) from Eq.(53) determines in terms of f 

where the function of integration (f){x) is independent of y. However, from 
Eq.(31), one may set 0(x) to zero using a gauge transformation with u!^{x^). 
Further, Eqs.(53) and (55) imply the same boundary conditions 

d,iB^h^-ln = fe^e^^(-«)Too(|/«) (57) 



3- 



y=ya 3-^5 
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Inserting Eq.(56) into Eq.(57) then gives 

o'f{x\y^)^t^Tooiya) (58) 

where in our static approximation V^. While is not gauge invariant 

in the bulk, we saw that it was gauge invariant on the branes, which is why 
its value on each brane is determined by the physical quantities Too(ya)- 

If one now inserts Eq.(56) back into Eqs.(53) and (55), one sees that these 
eqTiations are identically satisfied and so Eqs.(53) and (56) have no further 
content. Thus rather than determining and separately, Eqs.(53) and 
(55) determine only the gauge invariant combination Eq.(56). 

To check the solution of Eq. (44) , we first note that d^Tij is second order 
and may be neglected. Using Eq.(54) and the coordinate conditions Eqs.(49) 
and (51), and the fact that 

A"h,-^6{y-ya)h, = (59) 

one sees that Eq.(44) reduces to 

- ■^Too{ya)Siy - ya) (60) 

The boundary conditions imphed by the right hand side of Eq.(60) are thus 
identical to Eq.(57), and inserting Eq.(56) one sees that Eq.(60) is identically 
satisfied. 

We will see in the following section that the remaining Einstein equations, 
Eq.(45) uniquely determines /i^-^ so that we have found solutions to all the 
Einstein equations. The undetermined function, /^(x',y) off the branes, is 
the remaining gauge freedom. However, note one cannot set h^^ = every- 
where (as is conventionally done in other analyses) as Eq.(50) would then 
imply /'^ is constant in y, which would be inconsistent with the boundary 
conditions Eq.(58) (which are in fact gauge invariant). We will see that 
Eq.(58) contributes a significant term to the static gravitational potential. 
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4 Solution For hj^ 

The remaining Einstein field equations, Eq.(45), can be obtained by taking 
the TT part of Eq.(54). We find 

C-dl - 2A'd, + ^e^^n^)/.5^ = E^r(^")^(^ - (61) 

and hence /i^^ obeys the boundary conditions 



2A 

= (-^)"^^n^«) ; «=1,2 (62) 



The static potential is obtained from hQo{x,ya) where 



The corresponding source is then 



which in the static limit reduces to 



= It^o (65) 



To solve Eq.(61) we Fourier analyse h 



TT 
ij 



h'J{x\yo)= d'pe^^^hjnP\y) (66) 



In the bulk then hj^{p, y) obeys 

{\dl - 2A'95 + \e'^m^)hj^{p, y) = (67) 

where rn? = —p^ = pl—j^. The solutions of Eq.(67) are Bessel and Neumann 
functions 

}^;'{p\ y) = e'f'y[Mp)MO + BMm)] (68) 
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where 



(69) 



and m/P is short hand for (m^//5^)^/^. The boundary conditfons Eq.(62) 

One finds on the branes 



determine A^j and Bij 



h^ip; yi) = 



3(3MI 



-Too{yi) H p; -?oo(y2)J 



where 



D = 



D 



iVi(6) iVi(6 



m 



^1(6) ^1(6) 
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AT- 



12 



■/2(6) 

. ^2(6) 
6^1(6) 



Mil) 

"A^2(6) 



iVi(6) " 

Niiii) 



MCi) Mil) 



(70) 

(71) 

(72) 
(73) 
(74) 

(75) 



and 



2e2/^^%iV2i(6,e2)^ . . , A^22(^i,6) ^ . 

-/oo(yi) H ^ -?oo(2/2)J 



3/3M| 



(76) 



where 



A^2i(ei,6) = iVi2(6,ei) ; iV22(ei,6) = Arn(6,ei) (77) 

It is useful to consider /iqq^(j9, i/q,) as a function of a complex variable z = 
. In looking at the analytic behavior in z, the logarithmic branch cuts 
in the Neumann functions cancel in the differences such as A^i(Ci)/<^i(Ci) ~ 
Ni{i2)/ Ji{i2)- Poles can arise from a number of sources. Thus in Nu, a pole 
might occur at the zeros of Ji(^i), but this is actually canceled by the zeros 
in Ji(^i) appearing in D. Similarly the zeros of Jii^'i) in ^ii are canceled. 
Thus the only poles that occur are the pole at = (arising e.g. from 
N2/ J2 ~ l/m^ in A^ii) and when D vanishes, i.e., at 



D{ml) = 



(78) 
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where = {iT^n/P'^y^'^- To find the residue at m^, we expand -D(m^) around 
the pole position 



dD{m'^ 



[m^ — ml) + . . . 



and differentiating the Bessel and Neumann functions in D one has 

dDim'^) 1 1 



dm? 



2(5mJl{i 



2,.^['^l(el)iv^(el)-iVl(el)^^(el) 



(79) 



(80) 



2(5mJl{i2) 

Using the two Bessel function Wronskian identities 



[Ji{^2)N[{i2)-N^{^^)J[{^^)] 



JiN[ - NiJ[ = — = N,J2 - J1N2 



Eq.(79) reduces to 



-]) (m'-m^) + ... 



(81) 



(82) 



To obtain the residue at the pole, we need also the numerator A^n evaluated 

at m^ = ml 



6 JKCi 



■[A^2(6)^i(6)-A^i(6)^2(6)] 



and using Eq.(81) this reduces to 



Hence the residue at m^ is simply 



nml Jl{ml/p^) 



2 iVn 



3pMi D 



4/3 



3M? 



5 1 



Jf{mn/P) 
jf((m„//3)e^J'2) 



(83) 



(84) 



(85) 



The residue at m^ = can be obtained by examining the limit when both ^1 
and ^2 approach zero in A^ii(^i, ^2)/-D. One finds 



Rn 



A(3 



3M| 1 - e-2/3j/2 



(86) 
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where for the Randall- Sundrum model one may neglect the e^^*^'^^ ~ 10"^^ 
in the denominator. A similar analysis to the above holds for the other three 
terms in Eqs.(70) and (76). 

One may now cast the results for h^Q^z^Ha) in a more convenient form. 
Using the asymptotic forms of Bessel and Neumann functions, one can see 
that h^Q{z,ya) falls hke 1/z for \z\ on a large circle. Hence integrating 

g{z) ^ M^ifi^ (87) 

over a large circle in the complex plane we can express h^Qim?, Uq) as a sum 
of poles with the residues calculated above. One has then 

= + e-^^^^roo(y2)] (88) 



and 



hl^{y2) = -^^^^ooiyi) + e-'Py^T,,{y,)] (89) 



Ml 



The first term of Eqs.(88) and (89) contributes to the Newtonian potential 
(since = — — > in the static limit) while the other term gives 5D 
corrections to the Newtonian theory. 

While Eq.(78) is a transcendental equation, one can obtain the positions 
of the poles analytically in certain limits. Thus if = nT-n/P ^ 1 but 
^n^^^^ ^ 1 (i.e. ^ 10^^^) then inserting in the Bessel function asymptotic 
forms in the ratios of Eq.(78) gives 

tan{iyy^-^-^)^--{^f (90) 

which can be solved by iteration to give 

^ = [(^ + \)^ + enfe-^^y- ; « 1 , U^'' » 1 (91) 
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where 

6„-^[(n+^)7re-^^f (92) 

(We have included the first order correction as in some expressions the 
leading term can cancel out) . The residues at the poles can then be calculated 
in this hmit. Thus for the Too(yi) term one finds for 



P ^ 2 7Vn(6,6) 



3PMi D 

the result 



(93) 



277- 

Rn = -^^ne-"'' (94) 



and the contribution to the scalar potential is 



f'oo = E ^^^^^Too{y^) (95) 

o n=l " 



Since the poles are very dense 

Am„ = m„+i - m„ = pne'^^^ (96) 

one can approximate Eq.(95) by converting the sum to an integral (m^ = 
—p^ — —p^ in the static limit) 



27r rf^ 
'^00 - r>n/r3 



I ^^^^ 

Jq P ~r m 



3M| Jo f + ml 

where we have cut off the integral at /? since = m„//5 < 1. Returning to 
coordinate space this yields 



where tiiq is the mass of Too(yi). For r ^ i.e. for distances large 

compared to the warping parameter 1/(3, this is a correction to the 
leading Newtonian potential. 

Eventually, for sufficiently large n, ^„ becomes large (i.e. n > 10-*^^), and 
the poles from Eq.(78) occur at 

7X71 (3 
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/3r 



Then one finds for tliis contribution 



'"J»^--5|E^^^««fe.) (100) 



or in tlie continuum approximation 



In coordinate space one finds 

and in the limit r -C 1//3 one finds a 1/r^ correction to the Newtonian 
potential 

^" M - 

Eqs.(98) and (103) agree with results obtained in (although there, a fine 
tuning of matter is needed on the second brane in order to get a consisitent 
solution of the Einstein equations). 

One may carry out a similar analysis of the other three terms in Eqs.(70) 
and (76) {N12/D, N21/D, and N22/D) and these results will be discussed 
further in the Appendix. In Eq.(97) it is conventional to extend the integral 
down to TTin = 0, and think of the continuum of poles as reaching down 
to = without a gap. Actually, as can be seen from Eq.(91), the first 
discrete pole occurs at 

m„ = {^n(3)e-^y' (104) 

The size of the gap depends on the model. Thus for Randall-Sundrum one 
has 

nil ~ {^7ip)e-^y' ^ (10^^G'eV)(10-^^) = ITeV (105) 

since 

P ^ Mpi (106) 

On the Planck brane ?/i = 0, a TeV of energy is negligible (since masses are 
of order Mpi). On the TeV brane 1/2 however, it is sometimes argued that 
one should not consider phenomena > ITeV. In this case one would neglect 
the Kaluza-Klein modes. 
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5 Newtonian Potential 



The static Newtonian potential is the 1/r terms of hoo{x\ya) of Eq.(63). 
These arise from the poles in momentum space at = 0. As discussed in 
Sec.4, these poles occur in from the fact that the numerator functions 
Nij go as Nij ~ as — > due to the N2{^a),o: = 1,2 terms, while 

the denominator function goes as D ~ 1/m^ leading to a net 1/m^ term 
for small m^. As seen from Eq.(58), f^{x\ya) is totally a 1/m^ term in 
momentum space. One can thus pick out the = pole contributions on 
the two branes 

- -^i^2iToo{y.) + e-'^y^TM] + ^i^^TM (107) 

and 

C(y2) = -^i^ATooiy.) + e-'^y^TM] - ^.^,TM (108) 

In Eqs.(107) and (108) the first bracket is from Hqq and the second is from 
f^. The stress tensor T^j arising from the matter Lagrangian jC^n is 

= -L^ (109) 

where for a point particle on brane i/a 

J^m^ = mo y dTu'u^gij{x\ya)5^{x' - x\t)) (110) 

and = dx^dr with dr^ = —gijdx^dx^ . For our metric, gij = e~^^g}j where 
in the linearized approximation g'ij — rjij + hij. Thus defining 

dx^ 

dr'^ = g'ijdx'^dx^ ; u' = — (111) 

(Jj I 

Cm reduces to 

C^ruc, = rhaiva) J diu'ii^ gij5^{x' - x\f)) (112) 

where 

fh{y) = e-^^y^rriQ (113) 
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showing the usual result that if mo is of Planck size, the effective mass seen 
on the TcV branc ?/2 will be of TeV size. The Lagrangian of Eq.(112) will 
then correctly give rise to the (linearized) geodesic equation governed by 

9ij ~ Vij ~l~ ^ij- 

Returning to Eq.(109), the stress tensor is 

r^' = -^mo / dTu'u^5\x' - x\t)) (114) 

and in the static approximation, 

vP'^e^ ; u'^Q (115) 

one has 

TOO ^ e^Vo5'(r - r{t)) (116) 

so that 

Too(y«) = e^^^y-^m{y^)5\r - r{t)) (117) 
The interaction potential between the two particles may be defined by 



V 



= -/dV.„.„, (118) 



where the total is 



Cm = / c?™V^e-2^(^-)(7y,,- + hij{y^))5\x - x^ir)) (119) 



Hence in the static limit 



- j dhCm = - E^o« / dT(^fe-'^^y"\7^oo + hM)5\x - xjr)) 

(120) 



Since vP^dr = dx^^ and 



dxl^ 1 _ e 



A 



121) 



one has 

- j dhCm = 5] '^a(-?7oo - hoof'' (122) 



Expanding to first order gives for tlie interaction potential energy 



^ = -^$^^"^oo(x°,2/a) (123) 

a 

Inserting Eq.(107) and (108) and returning to coordinate space (m^ = —jf) 
one gets for tlie Planck brane the contribution 

P 1 _ _ 4 _ _ 
" "Mllr^^^^'^ ^ 3"^i"^2] (124) 

where rhi is the mass of a second particle on the Planck brane, m2 a mass 
of a particle on the TeV brane (m2 = e~^^^m2o), and r is the 3D distance 
between the particles. Note that the fact that m2 is separated by additional 
distance in the fifth dimension {1/2 — yi = vrp) does not enter in r. 

We see from Eq.(124) that if the two particles are on the Planck brane, 
Eq.(126) correctly reproduces the Newtonian force law with 

(the conventional value for the Newton constant in Randall-Sundrum the- 
ory). The f'^ contribution correctly changes the 4/3 factor in the first term 
of Eq.(107) to 1. However, if one particle is on the TeV brane, the Newton 
constant is modified by an extra factor of 4/3, since the /"^ factor does not 
contribute. (The fact that matter on the TeV brane changes gravitational ef- 
fects seen on the Planck brane has previously been noted in |2I] in a different 
connection.) 

For the potential energy seen on the TeV brane we use Eq.(108) in 
Eq.(123). One finds now from Eq.(117) that 



3M| Svrr 



■msm'ae^'^^^ (126) 



where m'g is a second particle on the ?/2 brane. The interaction energy between 
m2 and rhi particles is as before as is the Newtonian potential between two 
particles on the TeV brane, 777-2 and 7772, arising from h^ff (aside from the 
peculiar 4/3 factor). However, the term gives an additional contribution 
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to V{y2) scaled by e^"^*^^^^ (the factor from Eq.(117)) which would produce 
an anomolously large additional contribution. (Recall e^'^^ ~ 10^^ in the 
Randall- Sundrum model to account for the gauge hierarchy problem!) Thus 
the theory does not appear to give sensible results on the TeV brane. 

6 Conclusions 

In this paper we have examined the gravitational forces between point parti- 
cles in the static limit in the two brane Randall-Sundrum model. In contrast 
to previous analyses, we have chosen gauge conditions (coordinate frames) to 
solve the field equations that maintain the / Z2 boundary conditions, and 
hence produce no brane bending effects, and we also examine forces between 
particles on both branes, not just the j/i = brane. A convenient technique 
for solving the field equations is to introduce for the 4D generalization of the 
ADM decomposition [2H1 1211 metric perturbation 

hij = hj^ + hjj + hi^j + hj^i (127) 

where d'H(^ = = d'h^, ri'^hj^ = 0, and r/^^Tig = ^ 0. The hf^ 
contain the Kaluza-Klein modes while both the hf/" and hj, contribute to 
the static Newtonian potential (with pole at p ^ = in momentum space). 
One finds that a particle on the yi = brane sees a Newtonian force from 
another particle on either the yi brane or the 1/2 = ^rp brane but with different 
Newtonian constants: Gn = P/SnM^ and Gn = P/QnM^ respectively. The 
difference arises from the fact that the component of /loo enters with 
opposite sign for yi and y2 particles, as seen in Eqs.(58) and (63). (The fact 
that matter at 7/2 effects matter at yi differently from other matter at yi was 
also noted in [21] in another connection.) Note that the contribution is 
precisely what is needed to give the conventional value Gn = P/8nM^ on 
the yi brane. 

A curious feature of the potential Eq.(124) is that the force depends only 
on the 3D distance, and is independent of any y separation. It would be 
interesting to see if this produces any causal questions i.e. if one jiggled the 
mass on 7/2, how long does it take for the effect to become noticeable at the 
yi particle, a question involving dynamical rather than static solutions. 

A more serious problem is the force seen by two particles on the TeV brane 
y2- One sees from Eq.(126) that there is an attractive term arising from the 

contribution which is (^(e^^^^) ~ (10^^)^ larger than the normal gravity 
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and this occurs after one has correctly rescaled the 1/2 masses to TeV size 
(as one normally does in the RS model). Thus one does not recover normal 
Newtonian gravitation in the static limit on the TeV brane. All analyses up 
to now have neglected the Goldberger-Wise scalar field. Including it in might 
in some way cancel out the anomolous e^^^^ factor in Eq.(126). The analysis 
including the scalar field is much more complicated than the calculation given 
here. Preliminarily we find equations for at the brane positions to be 
given in terms of the scalar field, 0, by 

^'fiVa) = -Mf^'Too - -JJfWM - m] (128) 

Where 6(f) is the pertubation of the scalar field about the vacuum value 0o- 
The perturbation 6(f) satisfies a complicated equation of motion given in terms 
of the bulk potential V{(f>o) and its derivatives with respect to (f> as well as 
the brane potential Va{(j)o)- 

e^^U^6(f) - AA'6(f)' + 6(f)" - 6(f)V"{(f)o) - V\(f)o)h^^ 

= Y.^{y- 2/«)(^K(</'o)/i55 + v::{(f)o)6(f)) (129) 

a 

which is not decoupled from The remaining equations for and ^155 are 

^55 



At present it is not obvious that once the solution of this system is found 
(and then substituted into Eq.(128)) that this will rid the system of the 
exponentially large Newtonian force felt on the brane at 7/2. 

We briefiy compare our analysis with some of the previous calculations for 
the static gravitational potential. In Ref.j2I] it is assumed that in Gaussian 
coordinates 

^m5 = o ; /i = 0,1,2,3,5 (132) 



3A 



-d f 



3 A' Ml' 



(130) 



-^yi<f>o)h, + ^V\(f)^)6(f) 



Ml 



(131) 
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there is no brane bending, and brane bending occurs only when one adds the 
coordinate conditions 

d'hij =0 = hi ; ij = 0,1,2,3 (133) 

One can easily check, however, that the extra condition h^r, = of Eq.(128) 
cannot be achieved without introducing brane bending. Thus to achieve 
^55 = 0, we see from Eqs.(28) and (50) one requires 

Ux\y) = -^f{x\y) + Mx) (134) 

where the function of integration 05 (x) is independent of y and f'^{x\y) is 
the value of in the frame of Eq.(20). On the branes, therefore f^iva) is 
given by Eq.(58), and one cannot choose 05 (x) to make ^5 vanish on both 
branes. Thus if we choose 05(x) = -f'^{x\yi)/6A' (so that ^5(1/1) = 0) then 
^5(^/2) is proportional to e'^^y^m{y2) [by Eqs.(58) and (117)] and so there is 
a huge amount of brane bending on the TeV brane. (Alternately, the choice 
05(x) = gives by Eq.(30) that f'^{x\y) = but with brane bending on 
both branes.) This would presumably greatly modify the geodesic motion 
of particles on the TeV brane. Ref.j2Z| carries out the analysis in the frame 
of Eq.(128) assuming there is no brane bending in that frame. They define 
the gravitational potential by the diagram of two point mass stress tensors 
connected by a free field gravitational propagator. (The components will 
vanish for free fields.) However, to get a consistent solution they find it 
necessary to fine tune the matter on the ?/2 brane. In contrast, the analysis 
given here is valid for arbitrary matter on the yi and 2/2 branes. Finally we 
note that none of the previous discussions have analysed gravitational forces 
involving particles on the ?/2 brane which is where difficulties arose. 

Note Added 

After the completion of this work, it was brought to our attention that there 
exists a previous work that examines the corrections to Newton's law in a 
coordinate frame where h^^ 7^ and the branes are kept unbent [SO] • However 
we disagree with their results for the Newtonian potential. 
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A Appendix 

In Sec. 4 we calculated the Kaluza-Klein (KK) corrections to h^Q on the yi = 
brane in the case where both particles reside on the yi = brane. In this 
Appendix we will calculate the other KK corrections for the two cases (i) 
^1 -C 1, ^2 ^ 1 and (ii) S> 1, ^2 2> 1. These can be most easily found 
in terms of what we have already shown for h^{yi) due to the presence of 
Too(yi). These results were 

1) = E §^T,o{l) ; 1) = -^-^rn^ (A-l) 



— m. 

n 



and 



1) = E ; Kil, 1) = -^^^n (A.2) 



3M| 



where we have denoted the contribution to h^Q on the i'th brane due to 
matter on the j'th brane by lioQ{i,j) and (i), (ii) represent the two limits on 
,^1 and ^2 stated above. 

We can convert these into coordinate space using 

where the last approximation is true in the static limit. We then take the 
continuum limit where 

Am„ ^ drrin = Tr/^e-^^^ (A.4) 

Thus for Eq.(A.l) 
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where we have taken the upper hmit of the integral to be f3 since 1- 
After performing the coordinate space integral we get 

Upon integration we find 

Thus in the hmit r » 1/P we have a correction 

- e^Jw^^""^^^ ; 6 « 1 ; 6 » 1 (A.8) 

Similarly from Eq.(A.2) we get 

hoo ''HlA) - ^^^,e-^'Too{l) (A.9) 
which in the limit r <^ 1 //3 becomes a 1 /r^ correction 



^oo"^'^(l,l)^^^^roo(l) (A.10) 

We consider now the other corrections arising from Eqs.(88) and(89). We 
have for the correction due to Too (2) on the yi brane 

h^{i,2) = p^§^A^{l,l) (A.11) 



For case (i) we have 



37r^ 



Prom Eq.(90) tan{^2 - t) = -^li'^^i) and so 
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Thus 

^^)VV- (A.14) 



_ , W .1/2 



Substituting this expression into Eq.(A.ll) and after performing the coordi- 
nate space integration we find 

^^^^ = 6^^r\lp I dmnml/^e--^T,o{2) (A.15) 

We can calculate the integral in the limit r ^ 1/13 which gives a 
correction 

For case (ii) we need 



e" 2 (-1)" (A.17) 



After substituting this into Eq.(A.ll) and performing the coordinate space 
integral we find 



Wy2 



hoo^'Hh2) = 3^^3 ^ ;^(-l)"e-''-roo(2) (A.18) 

5 n 

Here since ^ 1 the sum is over m„ > /3 and since m„ = mrPe'^^'^ we 
require 

n > = AT > 1 (A.19) 



TT 

Thus Eq.(A.18) becomes 

3/3 00 

/^oo'^^^^(l,2) = 5](-l)'^e— ^--"^^700(2) (A.20) 

5 n=N 

Let 

n = Ar + m ; m = 0,1,2,... (A.21) 

Since 

Nnpre-''^^^y^ = /3r (A.22) 
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one has 



3/3i<2 



a - -fir °° 

C('^)(l,2) - V E(-l)"^""^'^"'"^oo(2) (A.23) 



j^(_l).^-..,.e-...^^^(2) = (A.24) 



5 m=0 

The sum is found to give 

Too(2) 

m=0 

Thus in the hmit r < 

C(^^)(l,2) = (-l)^^^^Too(2) (A.25) 
For the KK corrections on the y2 brane we have 

/ioo(2, 1) = e^^y-hoo{l. 2)^4^ (A.26) 
Thus from Eqs.(A.16) and (A.25) we find 

''^2' = ^1/2,5/2 ^00(1) ;/3r » 1 (A.27) 

/.S'(^^)(2,l) = (-l)^^^^^Too(l) ;/?r « 1 (A.28) 



Similarly the corrections for both particles on the y2 brane give 

^f(6)Too(2) 
'^f(e2)Too(l) 



/.-(2, 2) = (A-29) 



For case (i) we need 

Jl{i2) 2(3 

Hence 



(A.30) 
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and after going to the continuum hmit and performing the integrations over 
coordinate space and the mass spectrum we find in the hmit (3r ^ 1 

For case (ii) we have 

(6) 6co.2(6 - f ) 
which in the hmit /3r <^ 1 gives 



e^y' (A.33) 



^» ''(2,2) = ^^^Too(2) (A.34) 
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